Armchair graphene nanoribbons in presence of line defect 
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The characteristics of energy band spectrum of armchair graphene nanoribbons in presence of hne 
defect are analyzed within a simple non-interacting tight-binding framework. In metallic nanorib- 
bons an energy gap may or may not appear in the band spectrum depending on the location of 
the defect line, while in semiconducting ribbons the gaps are customized, yielding the potential 
applicabilities of graphene nanoribbons in nanoscale electronic devices. With a more general model, 
we also investigate two-terminal electron transport using Green's function formalism. 
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I. INTRODUCTION 



The fabrication of single layer graphene by Novoselov 
et alX in 2004 has opened a new era in the research of low- 
dimensional nanostructured materials. Graphene, the 
carbon allotrope with planar honeycomb lattice struc- 
ture, is a promising candidate of nano-electronic com- 
ponents owing to its exceptional electronic, thermal 
and transport properties^. It has been predicted that 
graphene sheet is a zero gap semiconductor—, but its 
behavior strongly depends on the boundary conditions 
when it is tailored into ribbon, flake or tube^. The 
sensitivity to the ribbon width, chirality, shape of the 
edges has allowed one to switch its semiconductor-like 
behavior from zero gaped to the finite one. Intensive 
researches have already been done on graphene nanorib- 
bons (GNRs) to explore the influence of edge topologj*^ 
on transport properties. Some attempts including chem- 
ical doping^-'™, application of uniaxial straiit^'^^, chemical 
edge modificationsii, incorporation of impuritj»i^, line 
defecti^"— are in focus of study on this system. But 
to realize the potential application of this material the 
control over transport properties needs to be clarified in 
a deeper way. 

To date, many theoretica l^^d^ as well as experimen- 
tal^)^ works have been done which reveal the fact that 
graphene nanoribbons (GNRs) with zigzag edges exhibit 
a metallic phase with localized states located on the 
edges, while armchair graphene nanoribbons (AGNRs) 
show metallic or semiconducting phase depending on the 
width of the ribbons22ri2i. This phenomenon is true only 
for clean ribbons i.e., without any deformation anywhere 
in the sample. But, the presence of impurity or defor- 
mation makes the system behave differently. In 2007 
Peeters et alt^ have shown that in presence of line impu- 
rity in graphene nanoribbon a gap opens up in the energy 
band spectrum. The system they considered was practi- 
cally coupled two graphene ribbons of different sizes sep- 
arated by a distance. Line defect also yields the possibil- 
ity of using graphene as a valley filter as demonstrated 
by Gunlycke et alr^. In a recent experiment topologi- 



cal line defect has been studied using scanning tunneling 
microscopei^ii^. Although the studies involving AGNRs 
have already generated a wealth of literature there is still 
need to look deeper into the problem to address several 
important issues those have not been well explored ear- 
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FIG. 1: (Color online). A schematic illustration of an arm- 
chair graphene nanoribbon, coupled to left and right leads, in 
presence of line defect (yellow circles) . The blue and magenta 
circles represent two different sub-lattices of the ribbon. 

lier, as for examples the understanding of the behavior 
of energy gap in metallic and semiconducting AGNRs in 
presence of line defect and also the dependence of this 
energy gap on the location of line defect. In the present 
work we mainly concentrate on these issues. Here we 
analyze the energy band structure of simple armchair 
graphene nanoribbons in presence of line defect within a 
tight-binding model. We show that depending on the lo- 
cation of a line defect in a metallic AGNR an energy gap 
may or may not appear in the band spectrum, while in a 
semiconducting AGNR the gap can be controlled. This 
phenomenon leads to the possibility of using AGNRs as 
nanoscale electronic devices. With a more specific model 
we also describe two-terminal electron transport using 
Green's function formalism to explore the results in real- 
istic cases. 

The structure of the paper is as follows. In Sec- 
tion II, we describe the specific model i.e., AGNR with 
side-attached leads and theoretical formulation to illus- 
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trate two-terminal electronic transport. The essential re- 
sults are presented in Section III which contains (a) en- 
ergy band structure of simple isolated armchair graphene 
nanoribbons in presence of line defect, and (b) transmis- 
sion probability as a function of injecting electron energy 
through the lead-AGNR-lead bridge system. Finally, in 
Section IV, wc summarize our main results and discuss 
their possible applications for further study. 



II. MODEL AND THEORETICAL 
FORMULATION 

Model : Figure [T] presents a schematic illustration of 
the model quantum system, where a finite size AGNR 
is coupled to two semi-infinite graphene leads with arm- 
chair edges. The blue and magenta circles correspond to 
two different sublattices of the ribbon, and, the yellow 
circles, representing the defect sites, are arranged in a 
line result a defect line. A unit cell of the AGNR is de- 
scribed by the dashed region which contains 2Ny atomic 
sites in our notation. 

Our analysis for the present work is based on non- 
interacting electron picture, and, within this framework, 
tight-binding (TB) model is extremely suitable for an- 
alyzing electron transport through such a two-terminal 
bridge system. The single particle Hamiltonian which 
captures the AGNR and side-attached leads gets the 
form: 

H = HagNR + -fflood + -fftun- (1) 

The first term i?AGNR denotes the Hamiltonian of the 
AGNR sandwiched between two graphene leads. Under 
nearest-neighbor hopping approximation, the TB Hamil- 
tonian of the AGNR reads, 

Hagnr = ^ Qcjci + ^ vii+i (clci+i + h.c.^ (2) 

where, e; is the on-site energy and is the nearest- 

neighbor hopping integral. The hopping integral 
is set equal to v or v' whether an electron hops between 
two ordered atomic sites or between two defect sites, and 
it is v" when the hopping of an electron takes place be- 
tween an ordered and a defect site, cj (c;) is the creation 
(annihilation) operator of an electron at the site I. 

The second and third terms of Eq. [T] correspond to the 
Hamiltonians for the semi-infinite graphene leads (left 
and right leads) and AGNR-to-lcad coupling. A similar 
kind of TB Hamiltonian (see Eq. [2]) is used to illustrate 
the leads where the Hamiltonian is parametrized by con- 
stant on-site potential eo and nearest-neighbor hopping 
integral vq. The AGNR is directly coupled to the leads 
by the parameters tl and th, where they (coupling pa- 
rameters) correspond to the coupling strengths between 
the edge sites of the ribbon and the left and right leads, 
respectively. 



Transmission probability: To obtain electronic 
transmission probability through the AGNR we use 
Green's function formalism. Within the regime of co- 
herent transport and in the absence of Coulomb interac- 
tion this technique is well applied. Using Fisher-Lee re- 
lation, two-terminal transmission probability T through 
the lead-AGNR-lead bridge system can be written aa^^, 

T = Tr [riGAGNRriiG'AGNR] (3) 

where, Tl and Ffl are the coupling matrices and G^^^j^ 
and G^gj^j, arc the retarded and advanced Green's func- 
tions of the AGNR, respectively. Now the single particle 
Green's function operator representing the entire system 
for an electron with energy E is defined as, 

G=[{E + nj)I-Hr^ (4) 

where, rj ^ 0^, H represents the Hamiltonian of the 
full system and / is the identity matrix. Following the 
matrix forms of H and G, the problem of finding G in 
the full Hilbert space of H can be mapped exactly to 
a Green's function Gagnr corresponding to an effective 
Hamiltonian in the reduced Hilbert space of the AGNR 
itself and we geli^^, 

Gagnr = [{E + if])I - /^agnr - Sl - ^nf^ ■ (5) 

Here, and are the contact self-energies introduced 
to incorporate the effect of coupling of the AGNR to the 
left and right leads, respectively. Below we represent 
explicitly how these self-energies are evaluated for the 
graphene leads attached to the AGNR. 

Evaluation of self-energy: In order to determine 
self-energies for these side-attached leads we follow the 
prescription addressed by Sancho et al^, where both 
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FIG. 2: (Color online). A schematic view of our model 
quantum system where AGNR and side-attached leads are 
sketched with discrete principal layers those are specified by 
integer numbers. 

leads and AGNR are sketched with discrete effective prin- 
cipal layers. These layers are defined as the smallest 
group of neighboring atomic planes and they allow only 
nearest-neighbor interaction between them. It effectively 
transforms the original system into a linear chain of prin- 
cipal layers^l as shown in Fig. [2l We label the principal 
layers in the right- lead as 0, 1, 2, . . . and in the left- lead 
as —1, —2, —3, ... and so on. The sample between the 
leads is denoted by C. Below, we describe elaborately the 
evaluation of the self-energy corresponding to the right 
lead, and, following this prescription wc also determine 
the self-energy for the left lead. Using Eq. 2] a set of 
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equations for the layer orbitals can be written as, 
[iE + iTj)I ~H 0,0] GoM = I + Ho,iGi,Q 

[(^ + z?/)7 - /fo.o] G2,o = /fJ,iGi,o + /7o,iG3.o 

[{E + ivi)! — HQfi]Gn,o = + Ho,iGn+i.o 

(6) 

where, we assume that .ffo.o = -f^i.i = -^2,2 = • • ■ and 
Hqj = Hi_2 = H2.3 = . . . . Here, Hij describes the 
Hamiltonian of ^-th principal layer, while Jf/^j+i corre- 
sponds to the coupling matrix between Z-th and (/-|-l)-th 
layers. The general expression of Gn,o reads, 



where, 



Gn.O — toGn-i^o + tQGn+1,0 

to = [iE + i7j)I -Ha^or'Hl i, 
to = [iE + i7j)I -Hofi]-'Ho,i. 



(7) 



(8) 



Substituting G„^o into the expressions of Gn-1,0 and 
Gn+i.o we can write, 

G„,o = tiG„_2,o + *iG„+2,o in>2) (9) 

and this process continues iteratively to repeat. After 
i-th iteration we have, 



G„,o = t.G„_2so + t'G„+2',o {n>T) (10) 



with. 



-1.2 



(11) 



The iteration is to be done until ti, ti < 5 with S arbi- 
trarily small. 

Using these t^'s we can determine the Green's function 
of a single layer in terms of the Green's function of the 
following or preceding one like, 

Gi,o-rGo,o and Go,o = TGi.o (12) 
where, T is the transfer matrix and it is defined as, 

T = {to + toil + totit2 + . . . + 40*1*2 ■ • ■ tn), 

f = (to + toil + totii2 + ...+ totlt2 . . . t„). (13) 

After some algebraic calculations we can write from 
Eq.H 



Go,o = [(E + iTj)I - Ho.o - Ho.iT]- 



(14) 



With this formalism, the surface Green's function of the 
left and right leads can be found as. 



9-1,-1 = [{E + iTj)I-Ho.o-Hl 



-2,-1^ 



where, .ffo.o and H-2,-1 are the Hamiltonians for a prin- 
cipal layer in both layer and the tunneling matrix be- 
tween two principal layers in the right-lead, respectively. 
The main advantage of this framework^^Ti^ rather than 
any other method is that here the number of iterations 
required for convergence is very small^i^. Finally, we 
get the expressions for the self-energies of the two leads 
as. 



— -ff-i,cfl'-i,-i-^-i,C' 



(16) 



where, ifo,i and Hi 2 are left lead-to- AGNR and AGNR- 
to-right lead coupling matrices, respectively. Using the 
above expressions of self-energies for the graphene leads 
we evaluate effective Green's function JTagnr with the 
help of Eq. [5] and then calculate two-terminal transmis- 
sion probability. 



III. RESULTS AND DISCUSSION 

In this section we present analytical results of energy 
band spectrum for isolated AGNRs and numerical re- 
sults computed for transmission probability through AG- 
NRs under conventional biased conditions. Throughout 
our analysis we set the on-site energies in the two side- 
attached graphene leads to zero, eg = 0, and in the 
AGNR 6/ = for ordered sites, while e/ = leV for de- 
fect sites. The nearest-neighbor coupling strength in the 
leads (wo) is fixed at leV, and the coupling parameters 
tl and tr are also set at leV. In AGNR, sandwiched be- 
tween two leads, we use three different hopping integrals, 
w, v' and w", and their values are fixed at leV, 0.7eV and 
O.leV, respectively. We fix the equilibrium Fermi energy 
Ep at zero and choose the units where c = e = h = \. 
The energy scale is measured in unit of v. 



A. AGNR without side-attached leads: Energy 
band structure and related issues 

To find the energy dispersion relation of a finite width 
AGNR we establish an effective difference equation anal- 
ogous to the case of an infinite one-dimensional chain. 
This can be done by proper choice of a unit cell (see 
Fig. [T]) from the nano-ribbon. With this configuration, 
the effective difference equation of the AGNR reads. 



{EI-S)ijj^T^,+i+Vi',-i 



(17) 



where. 



9o,o 



[(S + ir;)/-i?o,o--ffo,iT]-i (15) 



(18) 
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In the above relation, S and T correspond to the site- 
energy and nearest-neighbor hopping matrices of the unit 
cell, respectively, and X is the identity matrix. The di- 
mension of these three matrices is (2A^j, x 2Ny). Since in 
the nano-ribbon translational invariance exists along the 




FIG. 3: (Color online). Energy band diagrams of armchair 
graphene nanoribbons in absence of any line defect, where (a) 
Ny = 18, (b) Ny = 19, and (c) Ny = 20. 



spectra the three typical numbers (18, 19 and 20) of Ny 
are chosen only to make Ny in the forms 3n. 3n -I- 1 and 
3n -|~ 2, respectively, since the energy band structures 
of AGNRs are highly sensitive to these typical ribbon 
widths^. Here we set n = 6. From the spectra it is ob- 
served that for the particular case where Ny = 3?i + 2, the 
lowest conduction band and the highest valence band co- 
incides with each other at = 0, resulting zero energy 
gap in the band spectrum (Fig. [3ljc)). This indicates 
metallic phase of the AGNR. However, for the other two 
cases {Ny = 3n and 3n-|-l), a finite gap in the band spec- 
trum is obtained at fc^; = representing the semiconduct- 
ing behavior. In these three spectra (Figs.[3l^a)-(c)) since 
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FIG. 4: (Color online). Energy gap Ai5, at kx = 0, as a 
function of ribbon width A'^y. (a) Ny = 3n and (b) A'^y = 
3n -I- 1, where n is an integer. 

x-direction, we can write Tpj in terms of the Bloch waves 
and then Eq. [17] gets the form. 



{EX -£)= Te" 



(19) 



where, A = 3a is the 
unit cells, a is the leuj 
zene like ring. Solving 
dispersion relation {E 
As illustrative exam 
ergy band diagrams of 
widths when they are 



spacing between two neighboring 
5th of each side of hexagonal ben- 
Eq. [12] we find the desired energy 
vs. kx) of the armchair ribbon, 
pies, in Fig. [3] we present the cn- 
AGNRs for three different ribbon 
free from any line defect. In this 



FIG. 5: (Color online). Energy band diagrams of armchair 
graphene nanoribbons in presence of line defect, (a) Ny = 18, 
N = 5; (b) Ny = 19, N = 4; (c) Ny = 20, iV, = 6 and (d) 
Ny = 20, Ni = 5. Ni describes the location of a defect line. 



Ny is finite, the wavevector along j/-dircction becomes 
quantized and for each value of ky we get a E-k^ curve 
which results distinct energy levels in the E-k^ diagram. 
For large enough A^j,, energy gaps between these energy 
levels decrease sharply, and therefore, quasi-continuous 
energy bands are formed. The energy gaps between the 
conduction and valence bands, at fc^; = 0, of AGNRs with 
Ny — 3n and Ny = 3n -I- 1 strongly depend on the value 
of n i.e., ribbon width. To reveal this fact in Fig. 2] we 
present the variation of energy gap AE as a function of 
Ny for two different cases. It shows that the energy gap 
sharply decreases with Ny when Ny becomes smaller, but 
it eventually saturates to a finite non-zero value, though 
it is too small, for large enough Ny. So, in short, we can 
emphasize that a metallic phase is observed for AGNRs 



when Ny ~ 3n + 2, while the semiconducting phase is 
visible for the AGNRs with A^j, = 3n and 3n -I- 1. 

The results described above for clean nanoribbons i.e., 
nanoribbons in absence of any line defect have already 
been established in the literature, but the central issue 
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of our present investigation - the interplay between the 
existence of a hne defect, the width of AGNRs and the 
location of line defect has not been well addressed earlier. 

To explore it, we present in Fig. [5] the energy band 
diagrams for some typical AGNRs in presence of a line 
defect, where the defect sites are described by the yellow 
circles as shown in Fig. [TJ The location of a line defect is 
described by the variable Ni and we assign Ni = 1 for the 
edge of an AGNR. In Figs. Elja) and (b) the results are 
shown for the ribbon widths Ny = 3n and 3n + 1, respec- 
tively, and both for these two cases a finite energy gap 
around fc^; = is obtained which reveals the semicon- 
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FIG. 6: (Color online). Energy gap Ai5, at kx = 0, as a 
function of impurity position Ni when (a) Ny = 33, (b) Ny = 
34, and (c) A^^ = 35. 

ducting nature. The situation is somewhat interesting 
when the width of the AGNR gets the form 3n + 2. The 
results are shown in Figs. [S{c) and (d) where we choose 
Ny = 20 (= 3 X 6 + 2) and locate the defect lines at 
6 and 5, respectively. In one case a sharp crossing be- 
tween the energy levels takes place at kx = 0, results a 
metallic phase, while for the other case a finite energy 
gap opens up for this typical value of revealing the 
semiconducting phase. Thus, a metallic AGNR (width 
Ny = 3n -f- 2) can exhibit a metallic or a semiconducting 
phase depending on the location of a impurity line in the 
ribbon. In a metallic AGNR two types of states, metallic 
and semiconducting, for carbon chains exist which result 
these two types of conducting phases depending on the 
location of line defect, while a semiconducting AGNR 
contains only semiconducting chains, and accordingly, it 
does not provide the metallic behavior in presence of a 
line defect. 

The energy gap Ai?, across k^ = 0, strongly depends 



on the structural details i.e., the location of defect line 
in the nanoribbon. As illustrative example, in Fig. [S] we 
show the variation of /S.E as a function of impurity po- 
sition Ni for three different ribbon widths. It shows an 
oscillating behavior with the position of the defect line. 
For the semiconducting AGNRs (Figs.lHlJa) and (b)) the 
energy gap never drops to zero, while for the metallic 
ribbon (Figs. [Bfc)) /S.E vanishes when Ni becomes equal 
to 3p, p being an integer. These phenomena promote a 
design concept based on the structural details as semi- 
conducting devices with variable energy band gaps. 



B. AGNR with side-attached leads: Two-terminal 
transmission probability and ADOS 

Keeping in mind a possible experimental realization of 
the system, we clamp a finite sized armchair graphene 
nanoribbon between two ideal semi-infinite graphene 
leads (the left and right leads) making a lead-AGNR- 
lead bridge (sec Fig.[T|). Below we present our numerical 
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FIG. 7: (Color online). Transmission probability (green color) 
and average density of states (red color) for some typical arm- 
chair graphene nanoribbons in absence of line defect, where 
(a), (b) and (c) correspond to A^'j, = 3n, 3?! -I- 1, and 3n -I- 2, 
respectively. Here we set n = Q and take twenty unit cells. 

results for average density of states (ADOS) and two- 
terminal transmission probability through finite sized 
AGNRs under conventional biased conditions. 

In Fig. [7] we show the variation of transmission prob- 
ability, T, together with the average density of states 
(ADOS), p, as a function of energy E for some typical fi- 
nite sized AGNRs in the absence of any line defect, where 
(a), (b) and (c) correspond to the results for Ny — 3n, 
3n -f 1 and 3n -|- 2, respectively. Here we set n = Q and 
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choose twenty unit cells. Our previous analytical argu- 
ments for the AGNRs are exactly corroborated in these 
diagrams. A finite energy gap in the transmission proba- 
bility associated with the energy gap in ADOS spectrum 
is obtained when iV^ becomes identical to 3n and Sti + 1, 
as shown in Figs. W[a) and (b). This behavior empha- 
sizes the semiconducting phase for these ribbon widths. 
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FIG. 8: (Color online). Transmission probability (green color) 
and average density of states (red color) for some typical arm- 
chair graphene nanoribbons in presence of single line defect. 
In (a) and (b) the results are shown for Ny — 3n and 3n + 1, 
respectively, considering Ni — 5, while in (c) and (d) the re- 
sults are given for Ny = 3n + 2 when the defect line is placed 
at 5-th {Ni = 5) and 6-rd (A'^; = 6) lines, respectively. Here 
we take twenty unit cells and fix n at 6. 

On the other hand a gap less spectrum is observed when 
Ny becomes equal to 3n + 2 (Fig. El^c)), which indicates 
the metallic phase of the AGNR. 

Finally, we describe the results shown in Fig. [SI where 
we present transmission probability and ADOS for some 
typical AGNRs in presence of single line defect. In (a) 



and (b) the results are presented for A^^ = 3n and 3n -I- 1 , 
respectively, and for both these two cases transmission 
function shows a finite energy gap across E = exhibit- 
ing the semiconducting nature, with reduced amplitude 
compared to the spectra given in Fig. [T] The struc- 
tural dependence on the conducting behavior in metal- 
lic AGNR {Ny = 3n -I- 2) in presence of line defect is 
clearly visible from the spectra given in Figs. [HIc) and 
(d), where the defect lines arc placed in the 4-th and 
3-th lines, respectively. For the first case, it provides 
the semiconducting behavior, while in the other case the 
metallic phase is obtained, which perfectly corroborate 
our previous analytical findings. 



IV. CONCLUSION 

To summarize, we have investigated in detail the char- 
acteristics of energy band spectrum of armchair graphene 
nanoribbons in presence of line defect within a simple 
non-interacting tight-binding framework. The essential 
results have been presented in two parts. In the first 
part we have presented analytical results of energy band 
spectrum for isolated AGNRs. From our analytical re- 
sults we have analyzed that depending on the location 
of a line defect a metallic AGNR can provide either a 
metallic or a semiconducting phase, while a semiconduct- 
ing AGNR provides only the semiconducting phase with 
variable band gap. In the second part, we have discussed 
numerical results for transmission probability together 
with ADOS, keeping in mind a possible experimental re- 
alization of the system. We have shown that our numer- 
ical results exactly corroborate the analytical findings. 
Though the results presented in this article are worked 
out at absolute zero temperature limit, the results should 
remain valid even at finite temperatures (~ 300 K) since 
the broadening of the energy levels of the AGNR due to 
its coupling with the metal leads is much higher than 
that of the thermal broadenin g^^i^"^" — . 
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